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ABSTRACT 

We have used a high spatial resolution Chandra observation to examine the core mass distribution of the un- 
usually regular cD cluster Abell 2029. This bright, nearby system is especially well-suited for analysis of its mass 
distribution under the assumption of hydrostatic equilibrium: it exhibits an undisturbed, symmetric X-ray mor- 
phology, and a single-phase intracluster medium (ICM). From the deprojected temperature and density profiles 
we estimate the total mass, and the contributions of the gas and dark matter (DM) components from < 3" to ~ 3' 
(< 4.4-260/i7o kpc, 0.001 -0.1r v i r ). The gas density profile is not adequately described by a single /3-model 
fit, due to an increase in the density at the center (r < kpc, < 12"), but it is well fitted by either a double 

/3-model, or a "cusped" /3-model. The temperature data increase as a function of radius and are well-fitted by a 
Bertschinger & Meiksin profile and approximately by a power-law T(r) oc r aj , with ctj = 0.27 ± 0.01 . 

Using the fitted profiles to obtain smooth functions of density and temperature, we employed the equation of 
hydrostatic equilibrium to compute the total enclosed mass as a function of radius. We report a total mass of 
9.15 ± 0.25 x 10 13 /iyo M Q within 260h^ kpc, using the chosen parameterization of gas density and temperature. 
The mass profile is remarkably well described down to 0.002r v ; r by the Navarro, Frenk, & White (NFW) profile, 
or a Hernquist profile, over 2 decades of radius and 3 decades of mass. For the NFW model, we measure a scale 
radius r s = 540 ± 90hjQ kpc (« 0.2r v j r ) and concentration parameter c = 4.4 ± 0.9. The mass profile is also well- 
approximated by a simple power-law fit (M(< r) oc r Qm ), with a m = 1 .81 ± 0.04 (corresponding to a logarithmic 
density profile slope of -1 . 19 ± 0.04). The density profile is too shallow to be fitted with the profile described by 
Moore et al.. The consistency with NFW down to < 0.01r v ; r is incompatible with the flattened core DM profiles 
predicted for self-interacting DM (e.g., Spergel & Steinhardt), and thus contrasts with the strong deviations from 
CDM predictions observed in the rotation curves of low surface brightness galaxies and dwarf galaxies. This 
suggests that while CDM simulations may adequately describe objects of cluster mass, they do not currently 
account properly for the formation and evolution of smaller halos. 

Assuming that the cD dominates the optical cluster light within its effective radius (R e = 52", 76hjQ kpc), we 
observe a total mass-to-light ratio M/Ly ~ 12M Q /L Q at r < 20h^Q kpc, rising rapidly to > IOOM0/L0 beyond 
2OO/770 kpc. The consistency with a single NFW mass component, and the large M/L suggest the cluster is DM- 
dominated down to very small radii (< 0.005r v i r ). We observe the ICM gas mass to rise from 3 ± 1% of the total 

mass in the center to 13.9 ±0.4% at the limit of our observations. This provides an upper limit to the current 

-l 12 

matter density of the Universe, fl m < 0.29 ± 0.03 h 70 ' . 

Subject headings: galaxies:clusters:individual (A2029) — dark matter — intergalactic medium — 
X-ray:galaxies:clusters — cosmological parameters 



1. INTRODUCTION 

The large mass-to-light ratios (M/L) of galaxy clusters 
(M/L B w 200-300/170 M o /L , e.g., Girardi et al. 2002) in- 
dicate that they contain large quantities of dark matter (DM 
hereafter). The nature and distribution of DM is the current 
subject of much theoretical work in cosmology, with detailed 
simulations yielding different expectations for the amount and 
distribution of DM in cluster cores (e.g., Navarro et al. 1997; 
Moore et al. 1999; Dave et al. 2001). X-ray observations of the 
density and temperature of the hot intracluster medium (ICM) 
gas probe the mass of a galaxy cluster, under the assumption 
of hydrostatic equilibrium . Such data potentially provide con- 
straints on cluster DM simulations, and thus test DM theory 
(e.g., Evrard et al. 1996; Arabadjis et al. 2002; Sand et al. 2002). 

Prior generations of X-ray satellites have provided a wealth 
of cluster observations, from which we have begun to build 
a picture of large-scale radial temperature variations (e.g., Ir- 
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win & Bregman 2001; Ettori et al. 2002a), as well as two di- 
mensional temperature maps of disturbed systems (e.g., Briel 
& Henry 1994; Markevitch et al. 1998, 1999; de Grandi & 
Molendi 1999; Johnstone et al. 2002). However, detailed tem- 
perature and density profiles at the smallest scales (r < 0.1r v i r ) 
exist for only a few systems, such as Virgo (Nulsen & Bohringer 
1995), which exhibit irregularities in their cores. While gravita- 
tional lensing studies provide a unique and important probe of 
DM in cluster cores (Dahle et al. 2002; Natarajan et al. 2002; 
Sand et al. 2002), they generally cannot be applied to nearby 
systems, and they may also be contaminated by other sources 
of mass along the line of sight. The advent of the Chandra 
and XMM-Newton satellites now allows us to measure the ICM 
properties with simultaneous spatial resolution comparable to 
optical studies, thus providing a completely independent mass 
estimator over the same spatial scales. The main criticism lev- 
elled at X-ray studies is the veracity of the hydrostatic equi- 
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librium assumption in the actual clusters under study. Several 
groups have now obtained X-ray constraints on the DM profiles 
of clusters of galaxies, which are apparently consistent with the 
CDM paradigm (e.g., David et al. 2001; Pratt & Arnaud 2002; 
Ettori et al. 2002b; Schmidt et al. 2001; Allen et al. 2001; Mat- 
sushita et al. 2002; Arabadjis et al. 2002). Such systems are 
either too distant for a detailed analysis at < 0. lr vir , or they con- 
tain morphological disturbances indicating possible departures 
from hydrostatic equilibrium, especially at r < lOO/z^o kpc. Al- 
though simulations suggest that the X-ray analysis of such clus- 
ters is generally unaffected (Tsai et al. 1994; Evrard et al. 1996; 
Mathiesen et al. 1999), some authors argue that mass estimates 
will be biased at small radii in such systems, which may reflect 
the majority of clusters (see especially Markevitch et al. 2002). 

Abell 2029 (A2029) is a nearby (z = 0.0767), well-studied 
cluster of galaxies which provides an excellent opportunity to 
probe the dark matter content of a massive object. It has a very 
high X-ray flux and luminosity, as well as a hot ICM, indicating 
a massive system. We have previously analyzed the Chandra 
observations of A2029 in Lewis, Stocke, & Buote (2002, Pa- 
per 1 hereafter), presenting the temperature and Fe abundance 
data. It exhibits a very regular optical and X-ray morphology, 
and is an excellent example of a relaxed system with no ev- 
idence of disturbances (e.g., shock fronts, filaments, or "cold 
fronts") present in other systems. In a morphological analy- 
sis of 59 clusters, Buote & Tsai (1996) found the global X-ray 
morphology of A2029 to be among the most regular in the sam- 
ple. Though it contains a wide angle tail (WAT) radio source, 
there are no coincident X-ray "holes", such as those found in 
the clusters Hydra A (Nulsen et al. 2002), and Perseus (Fabian 
et al. 2002). There is no optical or X-ray spectroscopic evi- 
dence for a cooling flow, though the X-ray temperature drops 
to 2-3keV in the central 5" (Jh^\ kpc). Thus, this system is 
almost uniquely well-suited for analysis of its mass distribution 
since the hydrostatic equilibrium assumption should apply with 
high accuracy. 

In the current paper we present estimates of the gas den- 
sity and temperature profiles (§2), the total mass and the DM 
density profile (§3), as well as the gas mass (§4). In §5 we 
make a comparison with the stellar mass profile. We discuss 
the implications of our analysis and present our conclusions 
in §6. Throughout this paper, we assume a cosmology of 
H = 70 /i 7 o km s" 1 Mpc" 1 , Q m = 0.3, and A = 0.7, implying 
a luminosity distance to A2029 of 347 Mpc and an angular 
scale of 1.45 kpc arcsec" 1 . 

2. OBSERVATIONS AND DATA REDUCTION 

A2029 was observed by the Chandra observatory for 20ks 
on the ACIS S3 chip at a focal plane temperature of -120 C (see 
Paper 1 for details). We have reanalyzed the data in a similar 
manner to Paper 1, but using the latest version of the Chandra 
calibration (CALDB 2.15). Briefly, the data were processed to 
mitigate the effects of charge transfer inefficiency (CTI) using 
the ACISCtiCorrector.1.37 software 3 (Townsley et al. 2000). 
We then subtracted the available source-free extragalactic sky 
background maps 4 using the make_acisbg software created by 
Maxim Markevitch (see, e.g. Markevitch et al. 2000). We have 
also applied the latest correction to the ARF files to account 
for an apparent time-dependent degradation in QE (using the 

3 Available from the Chandra contributed software page at http://asc.harvard.edu/cont-soft/soft-exchange.html 

4 The background maps have been CTI-corrected in exactly the same manner as the source data. 

5 See http://cxc.harvard.edu/cal/Links/Acis/acis/Cal_prods/qeDeg/ index.html for a description of the effect and links to the software. 



"corrarf" routine 5 ). This correction has a small but significant 
effect on our fitted temperatures at larger radii, resulting in sys- 
tematically lower values than our previous analysis (see §3.4 
regarding the effect on the mass profile). 

2.1. Binning and Spectral Analysis 

We extracted spectra in concentric, circular annuli centered 
on the peak of the X-ray emission. We have used a different set 
of annuli than that used in Paper 1, to optimize the constraints 
on the shape of the mass profile in the center as well as the slope 
of the temperature profile at all radii. This results in a total of 7 
annuli, though this choice does not affect our results (see §3.4). 

Our primary analysis relies on spectra in the energy range 
0.7-8.0 keV, but we have explored the effects of using a lower 
limit of 0.5 or 1.0 keV (see §3.4). Using XSPEC, we fit the ex- 
tracted spectra with the APEC plasma emission model absorbed 
by Galactic neutral hydrogen. We adopt the weighted-average 
Galactic value of N H = 3.14 x 10 20 cm obtained from the 
W3N# HEASARC tool. The model normalization (from which 
we calculate the electron density, n e , and gas mass density, p g ), 
the gas temperature, T g , and the Fe abundance were allowed to 
be free parameters in the fits, with all other elements tied to Fe 
in their solar ratios. To properly recover the three-dimensional 
properties of the X-ray emitting ICM, we have performed a 
spectral deprojection analysis using the XDEPROJ code of Buote 
(2000). We start at the outside working our way to the core in 
an "onion-peeling" method which accounts for the cumulative 
projected emission from the outer annuli. For details of our de- 
projection technique see Buote (2000) and Buote et al. (2002b). 
To estimate the uncertainties on the fitted parameters we simu- 
lated spectra for each annulus using the best-fitting models and 
fit the simulated spectra in exactly the same manner as the ac- 
tual data. From 100 Monte Carlo simulations we compute the 
standard deviation for each free parameter which we quote as 
the "la" error. 

2.2. Azimuthally Averaged Gas Density and Temperature 

In Figure 1 we present the radial gas density and tempera- 
ture data for A2029 (left and right panels, respectively). The 
horizontal bars indicate the widths of the annuli used, and are 
not error bars. Note the very small statistical errors in both 
p g and T g , which indicate the precise spectral constraints ob- 
tained from these data. Since we wish to fit smooth, parame- 
terized functions to the data (see below), we choose to assign 
an emission-weighted effective radius, r, to each annulus ;' (see, 
e.g., McLaughlin 1999): 



r, = [(r^W^ (1) 
While r should strictly be iteratively calculated from the fit- 
ted density profile, it is easy to show that eq. 1 is an excellent 
approximation to the correct r over the entire profile (e.g., Ettori 
etal. 2002a, see also §3.4). 

3. TOTAL ENCLOSED MASS ESTIMATE 

We make the assumptions of hydrostatic equilibrium and 
spherical symmetry, such that the enclosed mass is 



M tot «r)=- 



_rkJg_ 
Gfim p 



dlnpg dlnTg 



dlnr d\nr 



(2) 



A2029 Dark Matter Profile 



3 



Radius [arcsec] Radius [arcsec] 

1 10 100 1 10 100 




1 10 100 1 10 100 

Radius [h'j kpc] Radius [h'j kpc] 



FIG . 1 . — Left Panel: Chandra radial gas density profile of A2029. For clarity, large open circles are centered on the data points (the smallest error bars are difficult 
to see in the logarithmic scaling). Horizontal bars indicate the sizes of the annuli used to extract spectra, and the limits of the spherical shells in our de-projected 
analysis. Overlaid are the best-fitting cusp model (solid curve), single /3-model (dotted curve), and double /3-model (dashed curve). Right Panel: Chandra radial 
temperature profile of A2029. Overlaid are the best-fitting Bertschinger & Meiksin model (solid curve), and simple power-law model (dashed curve). 



where k is Boltzmann's constant, G is the constant of gravita- 
tion, fi is the mean atomic weight of the gas (taken to be 0.62), 
and m p is the atomic mass unit. To obtain the instantaneous log- 
arithmic derivatives necessary to evaluate eq. 2, we fit parame- 
terized models to both the density and temperature. By parame- 
terizing the p g and T g data we derive a mass profile that may be 
smoother than the true mass distribution. This approach, there- 
fore, is best suited for interpreting average properties of M( < r), 
such as its radial slope and comparison with DM simulations 
(also smooth), which are the focus of the present paper. Key 
advantages of this method are that it is simple to implement, 
and the mass profile is straightforward to interpret in terms of 
the input p g and T g profiles. 



3.1. Temperature and Density Profiles 

We initially fit the gas density data with the ubiquitous /?- 
model: 

p g (r) = p g0 [l + (r/r c ) 2 ]-W\ (3) 

where p g0 is the central gas density, r c is the core radius, and 
-3(3 is the slope of the profile at r r c . The result is over- 
laid on the data as a dotted curve (Fig. 1, left panel). Due to 
a peak in the profile at < 17/^ kpc (the first 3 data points), 
the /3-model does not provide an acceptable fit (see Table 3.1 
below). 



Table 3.1: Gas Density and Temperature Fits 



Ps-Model (x'/dof) (3 r c \"} a p 

cusp 6.6/3 0.54 ±0.01 53.4 ±4.4 0.55 ±0.03 

1- /3 101.8/4 0.48±0.01 26.4±1.1 

2- /3 2.0/1 0.34 ±0.37 3.7 ±4.0 
0.76±0.14 53.2±9.5 ■■■ 

r r Model (x 2 /dof) TootkeV] r c ["] a T 

B&M 14.3/4 11.1 ±1.6 122.1 ± 125.3 0.36±0.05 

Power 19.8/5 ••• ••• 0.27 ±0.01 

NOTE - For the cusp model, we find p Sc = 6.6 ± 0.8 x 10~ 26 g cnf 3 . 
For the double-/3 model, p g i = 3.0 ±0.5 x 10~ 25 g cnT 3 and p g 2 = 
5.6±1.8x 10~ 26 g cnT 3 . 

We explored two additional models: (1) the 'cusp' model, 
which is a modified (3 model given by 

p g (r) = p g 2 3 ^ 2 - a ^ 2 (r/r c r^[l + (r/r c ) 2 ]-^ 2+a ^ 2 , (4) 

where p gc = p g (r c ), and the a p parameter allows a steepening 
of the profile at r < r c , and (2) a double-/? model (e.g., Xu et al. 
1998; Mohr et al. 1999) given by 

Pg( r ) = yftii+Pp, ( 5 ) 

where p g \ and p g 2 are each given by eq. 3. 

The double-/? and cusp models both provide satisfactory fits 
to the data (dashed and solid curves, respectively, Figure 1, 
left panel), though the reduced \ 2 is slightly improved for the 
double-/? model. We present the results of the gas density fits 
in Table 3.1. It is apparent that both the cusp and the double-/? 
models are sensitive to a break at w 53", and that the parameters 
for the first component of the double-/? model are not well con- 
strained. We have chosen the cusp model as our "reference" fit 
for the rest of our analysis for two reasons: ( 1) it provides a sim- 
ilar quality fit with two fewer free parameters than the double-/? 
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FIG. 2. — Left Panel: Total enclosed cluster mass, obtained from the B&M fit (open circles) and the power-law fit (open squares) to the temperature data. The 
cusp model for p g was used in both cases. Power-law fits to the mass points are overlaid on both data sets (solid line: B&M T g model, dashed line: power-law T g 
model). We have used large open symbols to identify the data points, as some of the error bars are barely visible in this logarithmic plot. The first data points are 
also enclosed with a large open diamond to emphasize the large additional systematic uncertainty at this radius (see §§3.2, 3.4). Right Panel: Total enclosed cluster 
mass (data points enclosed with open circles), overlaid with three different mass models: NFW (solid curve), power-law (dashed line), and M99 (dotted curve). The 
total enclosed gas mass is plotted as data points enclosed with open triangles. We have also overlaid an estimate of the stellar mass (dot-dashed curves, see §5.1). 
The lower curve assumes aM»/Ly of 1, the upper curve 12. 



model; (2) it yields a positive mass in the inner shell, which the 
(3 models do not (see §3.4). 

The temperature data are generally increasing, and are fairly 
well-fitted by a simple power-law T(r) =Ar aT , with otj = 0.27 ± 
0.01 (overlaid as a dashed line). The power-law fit misses the 
central data point entirely, and predicts a steeply rising slope at 
the outer point. A somewhat better fit to the data is obtained 
with a Bertschinger & Meiksin (1986, B&M hereafter) profile 
(which adds an additional free parameter): 

nr) = T 00 [r/(r+r c )r T , (6) 
where is the asymptotic temperature at large radii. We have 
overlaid the two fitted models to the temperature data in Fig- 
ure 1 (right panel). The B&M model (solid curve) provides 
an improvement in the fit to the data vs. the power-law model 
(dashed line) at the smallest and largest radii. The parameters 
of the fits are given in Table 3.1. From the table we see that 
the core radius is not well constrained, and that the temperature 
asymptotes to 1 1 . 1 ± 1 .6keV at large radii. Previous measure- 
ments of the temperature of A2029 at large radii (Molendi & De 
Grandi 1999; White 2000; Irwin & Bregman 2001), indicate a 
gas temperature of T g ~ 5 - 8 keV, suggesting a turnover in T g 
between 300-500 h^ kpc. Future work incorporating data at 
larger radius will thus require a more sophisticated temperature 
model. 

We note that as we found in Paper 1, the ICM of A2029 is 
apparently a single-phase gas at all radii. Even within a radius 
of 3", a single APEC model provides an excellent fit to the data, 
and there is simply no spectral evidence to support an additional 
component such as a second temperature, a cooling flow model, 
or a power-law. The B&M model provides a very good approx- 



imation to the T g data except for underestimation of the fourth 
data point. We cannot determine if this is a statistical fluctua- 
tion, a systematic error in our analysis, or real structure in the 
T g profile. The last case would introduce a correction factor to 
the mass in that region, otherwise the temperature, and thus the 
mass profile appears quite smooth (see §3.4). 

3.2. Mass Data and Fitted Profile 

At each radius i, we have calculated the total enclosed grav- 
itating mass M(< rf) according to eq. 2, using various density 
and temperature model fits. The logarithmic derivatives of p g 
and Tg are evaluated at F,, and T g (7j) is interpolated from the 
fitted Tg model. We estimate statistical errors gm, on the mass 
data as follows: For each Monte Carlo simulation j of the p g 
and Tg data (§2.1), we obtain a pair of fitted profiles (p/a, and 
Tfj tj ) from which we calculate a set of mass values M(< r,)j. 
We thus obtain 100 mass values at each radius, from which we 
calculate the standard deviation and report it as the "la" error. 

For the reference pair of gas models (the cusp model for p g 
and the B&M model for 7^), we present the total enclosed mass 
of A2029 in the left panel of Figure 2 (open circles). The profile 
has a nearly constant slope, with a slight flattening in the inner 
3 data points. At our final data point, we obtain a total enclosed 
mass of 9. 15 ±0.25 x 10 13 /i^ M within 260/1^ kpc. The mass 
calculated using the power-law T g model is also shown (open 
squares). We note that there is significant uncertainty in the 
mass at the innermost data point: while the temperature error is 
relatively high (~ 10%), the slopes of both p g and T g are even 
less well-constrained. In fact, using the double-/? model p g fit 
obtains a flat slope at this point, resulting in a negative mass 
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value (see §3.4). 

To analyze the shape of the mass profile, we fit parameter- 
ized models to the best fit mass values. To estimate errors on 
the parameters of these mass models, we obtain a mass fit Mf itj 
for each simulated mass data set M(< r)j. As above, we there- 
fore obtain 100 values of each parameter in a given mass model, 
from which we calculate a standard deviation. We have over- 
laid in Figure 2 (Left Panel) power-law fits to both mass pro- 
files. The reference mass data are well-fitted by a power-law 
M(< r) oc r Q, \ with slope a m = 1.81 ±0.04 (solid line), over 
the entire mass range. The mass data obtained from using the 
power-law T g model are slightly higher below r = \lh^\ kpc, 
and are also well-fitted by a power-law with a m = 1 .76 ± 0.03 
(dashed line). 

Although the power-law is a good visual fit to the mass data, 
we have examined whether \ 2 is reduced by instead using a 
broken power-law (BPL) model: 



M(< r) : 



A(r/r h ) a ", r<r b 
A(r/r b ) a °"', r>r h 



(7) 



where r b is the radius of the break between inner and outer 
slopes. If all the data points are included in the fit, then we 
obtain r/, = 5- 14 kpc, with a in = 1.1-1.5, and a ou , con- 
sistent with a, n given above (these results assumed a cusp p g 
model, and either the B&M or power-law T g model). However, 
if the inner data point is excluded, then the BPL fits are consis- 
tent with the single power-law fits given above. Since the inner 
data point is subject to large systematic error due to the chosen 
p g parameterization (as above, see also §3.4), we can conclude 
there is no evidence for a significant break in the logarithmic 
radial mass profile. 

To account for more gradual deviations from a power law, we 
have also fitted the mass data with the Navarro, Frenk, & White 
(1997, NFW hereafter) profile, p(r) oc [(r/r s )(\ + (r/r,) 2 )]" 1 , 
the Hernquist (1990) model, p(r) oc [r(l + r) 3 )]" 1 , and the 
Moore et al. (1999, M99 hereafter) model, p(r) oc [(r/r s ) L5 (l + 
(r/Vs) 1 - 5 )] -1 . Integrating these density profiles obtains mass 
models for fitting; analytic forms can be found in the litera- 
ture (see, e.g., Klypin et al. 2001; Loewenstein & Mushotzky 
2002). These fits are presented in the next section. 

3.3. The Dark Matter Distribution 

Assuming the cluster to be dominated by DM (a point we 
argue below), the fitted total mass profile corresponds to an im- 
plicit DM density distribution. For any power-law fit, puu oc 
r QM , where cxdm = ct m -3. For the power-law fit to the refer- 
ence mass profile we therefore observe a dark matter density 
slope ofa DM = -1.19 ±0.04. 

In Figure 2, (right panel) we again show the total enclosed 
mass profile (open circles). To compare with other clusters, 
as well as theoretical expectations, it is convenient to scale the 
radius in terms of the virial radius, r v ; r 6 . We calculate one pop- 
ular predicted value of r v ; r as a function of emission-weighted 
global temperature, using the form given by Neumann & Ar- 
naud (1999), eq. 9, normalized to the 2.0-9.5 keV band Mass- 
Temperature relation given by Mathiesen & Evrard (2001), 
which we convert to our cosmology. If we extract one spec- 
trum in this band from the entire region within 186" which we 
observe with Chandra, we measure a single-temperature fit of 

6 The virial radius is taken to be the radius at which the matter density is 200 times 

7 Our emission-weighted temperature measurement is in good agreement with the 
(2002) 



7.54±0.15 keV 7 . We thus obtain r vir = 2.71 ±0.42/*^ Mpc for 
A2029. We have plotted the upper axis of Fig. 2 (right panel) in 
units of r v ; r , which shows that we are examining the dark matter 
profile on a scale from < 0.001 -0.1 r vir . 

We have overlaid fits to the mass profile of A2029 from three 
different mass models: a power-law (dashed line), an NFW 
mass model (solid curve), and an M99 model (dotted curve). 
While the power-law model provides a good overall fit (x 2 /dof 
= 24. 1 /5), the NFW model is preferred (x 2 /dof = 1 1.8/5). The 
data are more closely approximated by the NFW profile except 
for a ~ 1.5cr difference at the innermost data point, which has 
additional systematic uncertainties (as noted above) rendering 
the discrepancy insignificant (see §3.4). The Hernquist pro- 
file fit (not shown) is nearly identical to the NFW fit (\ 2 /dof 
= 11.4/5). The M99 model does not provide an acceptable fit 
to the overall profile (x 2 /dof = 250. 1 /5), though its small radius 
slope is compatible with the inner 3 data points. However, if fit 
solely to these points, it falls well below the remainder of the 
mass profile, which is better constrained. We note that although 
the NFW and Hernquist models improve \ 2 significantly, the 
relative differences between these models vs. the power-law 
mass model are small (< 10% between 17 and 260 h~^\ kpc). 
Unlike the BPL model, which presents a sharp break, the grad- 
ual change in the logarithmic slope of the NFW or Hernquist 
profiles better quantifies the small deviations of the mass data 
from a pure power-law. 

For the NFW profile, we find a scale radius r s = 540 ± 
90h^Q kpc, and a concentration parameter c = 4.4 ± 0.9. This 
allows us to calculate the value of r v ; r expected from the NFW 
model (r v ; r = cr s ), for which we obtain r v i r (NFW)= 2.39 ± 
0.62/i7o Mpc, in good agreement with r v ; r as predicted by the 
Mass-Temperature relation given above. 

Several other authors find that CDM halos (e.g., NFW or 
M99 profiles) are also consistent with their X-ray cluster ob- 
servations at r > 0.1r vir (e.g., Ettori et al. 2002b; Pratt & Ar- 
naud 2002; Schmidt et al. 2001; Allen et al. 2001; Arabadjis 
et al. 2002), while gravitational lensing studies report conflict- 
ing results in cluster cores (i.e., Sand et al. 2002; Natarajan et al. 
2002). The Hydra A cluster exhibits p DM oc r~ 13 over a similar 
range of radii to our analysis (David et al. 2001), however, there 
are prominent interactions between the X-ray gas and the radio 
source in that system, implying significant deviations from hy- 
drostatic equilibrium. Thus the DM profile inferred for A2029 
(which is comparatively undisturbed in this regime), provides 
evidence that the Hydra A results are robust, and is an impor- 
tant independent confirmation of CDM predictions. 

3.4. Possible Systematic Errors in the Mass Profile 

We have explored several alternative data reduction and anal- 
ysis choices to investigate systematic errors in the mass profile. 
We have calculated density, temperature and mass profiles in 
four different sets of annuli: (1) the original 11 annuli used 
in Paper 1, (2) 13 annuli, obtained by dividing the inner 2 an- 
nuli of Paper 1 into 4 smaller annuli, (3) 5 annuli, using much 
larger bins, (4) 7 annuli, using a combination of the very small 
annuli from (2) in the center with the larger, higher S/N an- 
nuli from (3). Binning choice (4) is the reference choice for 
this paper, which shows the most detail in the core, and pro- 
vides the highest S/N for the temperature estimates at larger 

the critical density required for closure of the Universe. 

BeppoSAX analyses of both Irwin & Bregman (2000) and De Grandi & Molendi 
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radii. For all 4 choices of binning, the resulting mass profiles 
are not perceptibly different from the reference results. We have 
also analyzed the data in energy ranges of 0.5-8.0, 0.7-8.0, 
and 1.0-8.0 keV, and repeated the entire experiment using a 
MEKAL plasma model rather than the APEC reference model; 
in all cases the slope of the fitted mass profile is within the lex 
errors of the reference fit. 

The effect on the mass profile of our choice of eq. 1 to calcu- 
late F can be tested by instead using 2 extreme cases for the es- 
timate of F: assumed density slopes of p g oc r° and oc r" 3 , which 
easily bracket all the instantaneous slopes of the observed den- 
sity profile (see also Ettori et al. 2002b). We find that the mass 
profile does not vary perceptibly with the choice of F. 

Our reference analysis is to use the same Period C back- 
ground files as in Paper 1. Variations in the X-ray flux of 
the extragalactic sky between our target and the background 
templates may render them inaccurate. We experimented with 
modifying the background flux by changing the effective expo- 
sure time of the templates to ±20% of its nominal value. While 
this can affect the estimated temperature in the outermost an- 
nulus, the overall temperature profile does not vary beyond the 
lcr errors, nor does the slope of the mass profile. We have cor- 
rected the data using the "corrarf" routine (§2), which results 
in lower T g values than we reported in Paper 1, by up to 16% in 
the outermost bin. This has the effect of slightly flattening the 
fitted T g profile, as well as the mass profile (the uncorrected pro- 
file has overall slope a m = 1.91 ±0.03). The calibration of the 
ACIS instrument below 1 keV is ongoing, but it is reportedly 
now accurate at the w 10% level (see footnote 5). 

In Figure 2 (Left Panel) we show the total mass values ob- 
tained from each of the two temperature profiles shown in Fig- 
ure 1 (in both cases, we use the cusp fit to p g ). The mass values 
all overlap in their \a error bars, as do the slopes of the mass 
profiles. We have also explored the choice of the model which 
is fit to p g . The overall power-law fits to the mass profiles de- 
rived from the cusp and double-/? profiles are indistinguishable. 
However, we note that the double-/? model results in a negative 
(and unphysical) mass value in the central bin (due to its large 
statistical uncertainty, it does not affect the overall profile fit). 
The single [3 model (which is an unacceptable fit to the p g data) 
also obtains a very flat slope for the inner 3 data points which 
results in negative mass values over that region. To obtain a 
positive mass, eq. 2 requires a negative sum for the logarith- 
mic derivatives: inspection of Figure 1 shows that the slopes of 
neither p g nor T g are well-constrained at the radius of the inner- 
most data point. The cusp and double-/? models diverge in this 
regime, but the data do not distinguish between them. If the gas 
was in fact isothermal at this radius our mass estimate would 
be positive for all 3 p g fits, but without much higher S/N data 
(at higher spatial resolution) we cannot infer such a state for the 
gas. We note that exclusion of the innermost data point has no 
perceptible affect on the mass profile fit. 

The use of parameterized functions for p g and T g has the ef- 
fect of smoothing those data, in turn resulting in a smoother 
mass profile (§3.1). We have also calculated the mass profile 
by using a simple "point-to-point" estimation of the p g and T g 
slopes where the logarithmic derivatives needed for eq. 2 are 
taken to be the slope between each pair of adjacent points. The 
mass values are significantly different only in the fourth and 
fifth bins, which lie above and below the reference data points, 
respectively. Nonetheless, the average slope of the mass pro- 
file is not altered beyond the lcr limits of the reference fit (as 



one would expect, since the reference p g and T g models pro- 
vide good fits to the data). Aside from the innermost data point, 
the slopes of p g and T g are well constrained regardless of the 
parameterization, and the overall mass profile is robust. 

We have assumed spherical symmetry for A2029, though its 
mass distribution exhibits an average ellipticity of 0.5 (Buote & 
Canizares 1996) within w 1 .6^70 Mpc inferred from the X-rays 
and optical light. However, we have measured the spherically 
averaged mass distribution, which is quite insensitive to the el- 
lipticity of the system (including the specific case of A2029, 
Buote & Canizares 1996; see also Evrard et al. 1996). This 
approach is also useful for comparison with simulations which 
also present a spherically averaged profile (i.e., NFW, M99). 

4. GAS MASS, BARYON FRACTION, AND O m 

We have calculated the total gas mass M gas (r) in each spher- 
ical shell directly and summed them to obtain an enclosed gas 
mass profile (i.e., we do not interpolate the p g values from the 
fitted profile as we did for T g in the total mass calculation; the 
results with either method are unchanged). In Figure 2 (Right 
Panel) we have overlaid the total enclosed gas mass profile 
(open diamonds). We see that the gas mass at the center of 
the cluster is < 3% of the total mass, rising to 13.9 ± 0.4% at 
our last measured data point. This may safely be regarded as 
a lower limit, as simulations and other measurements suggest 
a higher asymptotic value for the gas mass fraction in clusters, 

W™) ~ 0.2-0.3/z~, /2 (e.g., Allen et al. 2002), while we only 
have observations at r < 0.1r v ; r . 

Assuming that A2029 contains proportions of dark matter 
and baryonic matter equal to their Universal values, one may 
place an upper limit on Q m with f gas . Given the Univer- 
sal baryon mass density from big bang nucleosynthesis cal- 
culations and recent observations of the deuterium abundance 
(Qb(BBN)), one finds fl m = fi B (BBN)/f B , where f B is the Uni- 
versal baryon mass fraction. Given n B h 2 (BBN) = 0.020 ± 
0.001 (see, e.g., Buries et al. 2001), h = 0.7 ±0.07, and f B > 

f S as > 0.139 ± 0.004h~l 1/2 from this work, we estimate Q m < 

0.29 ± 0.03/z 70 ' , in agreement with current estimates (for a re- 
cent review, see Turner 2002). 

5. MASS BUDGET OF THE CORE 

Massive elliptical galaxies such as NGC 720 (Buote et al. 
2002a), and NGC 4636 (Loewenstein & Mushotzky 2002), are 
apparently DM-dominated all the way into their cores. This sit- 
uation is likely to exist in a cD cluster such as A2029, except at 
possibly the smallest scales (rwl - lOkpc), where the cD stellar 
mass component may be more important (Dubinski 1998). We 
have found that the gas in A2029 is only a few percent of the to- 
tal mass in the core of the system. But are we observing the stel- 
lar mass (e.g., the cD)? Dubinski suggests the turnover between 
stellar and DM-domination should occur at a few tens of kpc; 
but although we are specifically observing this regime, we see 
no significant break in the mass profile (see §3.2). The outlying 
fourth Tg point and the central p g peak over a single /3-model are 
intriguing, but the balance of the data suggest that any transi- 
tion from stellar-dominated to DM-dominated regimes is quite 
smooth in this system. 

5.1. Stellar Mass 

We have estimated the contribution of stellar matter to the 
system, by using optical observations of the cD galaxy in 
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A2029 and its diffuse halo. Uson et al. (1991) conducted a de- 
tailed optical analysis of A2029, and observed that the surface 
brightness of the cD+halo was well-fitted by a de Vaucouleurs 
(1948) profile out to « 300". They report a total fl-band lu- 
minosity of 5 x 10 n /rj"ooL© for the cD+halo, which is 41% of 
the total light of the cluster within 260" radius. We assume the 
cD+halo to be the vast majority of the light within its effective 
radius R e = 52" (76fr^ kpc). Using the Hernquist (1990) model 
as an approximate deprojection of the R X I A law, we normalize to 
the total light of Uson et al. (1991), converting to our cosmol- 
ogy, and the V band, and obtain a luminosity profile. In Figure 
3, we present the total cluster mass-to-light ratio (M tot /L v ) pro- 
file of the core of A2029, based on the cD+halo light profile. 
It is consistent with a constant M tot /L v w 12M Q /L Q within 
20/170 kpc, and rises rapidly at larger radius. 

Dynamical stellar analyses estimate the total masses of galac- 
tic systems, but cannot directly measure the stellar mass. At- 
tempts to resolve this with population synthesis techniques 
yield highly uncertain results for the stellar mass-to-light ratio, 
M*/L v = 1 - 12M /L Q (Pickles 1985). Using these bounding 
values for M*/Ly we plot in Figure 2 (right panel) associated 
stellar mass profiles {dot-dashed curves) for the cD+halo sys- 
tem. The data indicate that at r < 0.005r vir (rs I5hjg kpc) the 
stellar mass can potentially dominate the system, depending on 
the highly uncertain assumed M*/L v . The stellar material out- 
weighs the gas mass up to a radius between w 0.005 -0.04r v ; r 
(« 15- lOO/z^Q kpc). If we assert a plausible value, M*/L v = 
5 M Q /L for the system, we would conclude that the cluster 
is DM dominated down to the smallest scales measured here, 
which is consistent with a single NFW mass component (and 
also with the lack of a significant break in the total mass pro- 
file, §3.2). Alternatively, the gas density peak over the /3-model 
at r < kpc, or the outlying T g data point may indicate the 
presence of a stellar mass component in excess of (or differing 
from) an NFW DM halo. 

5.2. Velocity Dispersion 

An interesting comparison can be made with the optical ve- 
locity dispersion profile of A2029. As noted by Dubinski 
(1998), IC1101 (the cD galaxy in A2029) is one of the only 
brightest cluster galaxies observed to have an increasing veloc- 
ity dispersion, ay, (Dressier 1979). We may regard the stars in 
the cD galaxy (from which ay is measured) to be simply an- 
other family of tracer particles (as the hot ICM gas particles) 
bound to the same gravitational potential. 

Thus, the effective velocity of the ICM gas particles may 
be expected to mirror that of the stars. We estimate the ve- 
locity dispersion for the ICM gas as oj CM = kT g / /im p . The 
temperature data yield <rf CM oc r °- 27±0 01 5 while the cD velocity 
dispersion profile of Dressier (1979), measured over the range 
2- IOO/zjq kpc reveals a v oc r 25 . This correspondence is not 
required for the condition of hydrostatic equilibrium, but the 
similarity of the two species is striking. A2029 is one of the 



few systems where the shape of both the X-ray and optical ve- 
locity dispersions can be measured in detail: the consistency 
between them is further evidence that we are observing a dy- 
namically relaxed system where all the mass components are in 
equilibrium with the gravitational potential. 

6. CONCLUSIONS 

We have analyzed high spatial resolution Chandra data of 
the A2029 cluster of galaxies, obtaining well-constrained ICM 
gas density and temperature profiles on scales of 0.001 -0. lr v ; r 
(« 3-260/170 kpc). Fitting smooth functions to these profiles, 
we obtain mass profiles for A2029, measuring the shape of the 
total mass profile at unprecedented accuracy to a very small 
fraction of the virial radius. Our results are insensitive to most 
details of the data reduction, but are closely tied to the well- 
constrained temperature and density distributions. 

We find that the shape of the inferred dark matter density at 
< 0. lr vir is consistent with the NFW parameterization of CDM 
halos, but apparently incompatible with that of M99, though we 
note that individual objects may be expected to show significant 
scatter from a mean DM halo profile (Bullock et al. 2001). 

The consistency of the NFW model with the mass profile 
all the way down to < 0.01r vir indicates that there is no need 
to modify the standard CDM paradigm to fit the DM distribu- 
tion in this cluster, consistent with X-ray observations of other 
clusters at larger radii (see §3.3). This result contrasts with the 
strong deviations from the CDM predictions observed in the ro- 
tation curves of low surface brightness galaxies (e.g., Swaters 
et al. 2000), and dwarf galaxies (e.g., Moore et al. 1999) which 
inspired the self-interacting DM model (Spergel & Steinhardt 
2000) to explain the relatively flat core density distributions 
in these galaxies. Hence in light of the good agreement with 
the NFW profile in clusters, the deviations observed on small 
galaxy scales do not seem to imply a fundamental problem with 
the general CDM paradigm. Instead, it is likely that the numeri- 
cal simulations do not currently account properly for the effects 
of feedback processes on the formation and evolution of small 
halos. 

Our analysis suggests that A2029 is dominated by a single 
mass (i.e., DM) component at all measured scales below 0.1r v ; r , 
or that any transition from a stellar mass dominated component 
in the cD and a DM component is quite gradual. We also ob- 
serve a rising gas fraction from < 3% to > 14% in A2029, ob- 
taining an upper limit to Q, m < 0.29 ± 0.03/z 7 q , consistent with 
other current studies. 

The authors wish to thank Amit Lakhanpal for an initial 
calculation of the mass profile. This work was supported by 
Chandra grant GOO- 102 IX. 
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